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The analysis o£ simultaneous laser data is investigated using 
the method of trilateration. Analysis of data rom 

SesenteS with selected baseline rates and statron positrons The 
use of simultaneous Etalon ^ST^tlateration 

hffal^t 9 ^“/regional defcrmation^studies with monthly 

LAGEOS American solutions between 3 1 


in the 1970s the precision of laser data was about 10 - 15 
Cen tions e was^oor he Sincf thaftime^he intentional laser ranging 

COTmunity W has expanded^ ^Precision^of the a 

eS^ot precision on LAGEOS below 2 centimeters and much of that 
data has a single shot precision below 1 centimeter. 

During the 1970s the mathematical methods to determine 

satellite positions using ^°^“^ e ^^ 0 ^ al Sr P rimn m advantage of 
term used for that method is f or bit It provides a means 

“i™Si^rion: s : rz «. 

lower other outside influence » ^uh *s satellite^, oSit'detemination 
pressure, and gravity models t trilateration include the need for 

methods. The raj « in separating rang, biases 

simultaneous ranging data, method does not allow for 

rar^fo^ 

ssss.r s-nr- — 

Our method of trilateration requires simultaneous data from 

four or more laser ranging stations. There ^^oess 

steps that must occur to prepare the data for analys • p lects 

7 v-p»-nr all of the polynomials have been determine , 

data tape. After all of t p Y“ second intervals from the 

simultaneous data points are created at 30 second 

^time inLrwals 
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within the selected passes that have fewer than 4 stations are deleted 
A summary output of the available stations and data yield is also 

seDarated** StSP ' . Finally ' the remaining LAGEOS tracking data is 

separated into two sections. One section is for Europe, North Africa 

Commonwealth of Independent States. The other section is for the 
Americas, Easter Island, and Huahine. 

ranaina geometricall y» each section of LAGEOS 

g g data is read in by the software. Apriori station positions and 
station velocities are used to compute each stations position at the 
middle of the month. The satellite position is determined for each 

to thTdat^ “ f? 91 ® dependent corrections are confuted and applied 
ia ai-a *■ rt — ?H ce al 1 the data has been corrected an iterative process 

satellite nofi> 6 ° bSerV ® d minus ^“"Puted values (O-Cs) from each 
satellite L v° n ^determined for each epoch along with an average 

criterion S °' CS ff* then * dited using a Siting 

9e Batelllte position is recomputed. A summary 

f sbatl ° ns mean and rms are written to the run summary file alSL 

ediL a ^ nUmber ° f observati °ns per station and the SLnber of 

edits. A sensitivity matrix is then generated. This is done by 

and da? 9 £ “! ta * 1 of ^ sets to each of the station position vector elements 

^fcha r ^ nin9 a effect each ° ffset has on the °- Cs for an of the 
pochs. This procedure is performed for each coordinate (ie 12 times 

“ t r* P ° S *, P 'then^inverted, 

ty.l t/- ! a partial derivatives of the O-Cs with respect to 

the adjusted station position. The O-Cs are then multiplied by the 

appUed^the^^at^ & ±S th ® n com P uted - These corrections are 

applied to the station positions and a new set of average satellite 

positions are computed and new O-Cs determined. Again these O-Cs are 

multiplied by the P matrix and reiterated. This iterative process is 

is P less t^ the d6lta rmS ° f tbe °' CS consecutive iterations 

one Percent, or six iterations have been reached. The 
total change of the station position is determined for both x y z and 

onli the'st^^ 11 ^ *"? hei9ht ’ The baseline lengths are determined for 
JS* Tt ® tatlo ” s inv °l ve d m the solution and saved in a database 
This database also contains the month and the number of simultaneous 
bservations between the two stations. A sample summary for a one month 
per od for a European solution is shown in Figure 1. A by-product of 

ifF WeT ^ th r- P ° int di8tribution a sample of whJch is shown 

1 ' Baseline rates are generated from the baselines and a 

bTwei^ are %i inear ° f 6aCh ° f the station combinations is computed 

Sat monti" 9 c “ ^ baseline with the number of observations^ring 

that month. Some sample baseline rates are illustrated in Figures 3 - 5 . 

t-o™ a i method of trilateration also lends itself to doing long 
a °^ tl0ns - 1” this method both the initial station position a2d 

under deve^ment 1 " 168 316 d ® termin ® d - method of analysis is still 

ThS req y irement for simultaneous ranging data from 4 stations 

hiS! ^ ■ r^ 9 raP ^fn COVerage over a broad area, which in turn requires 
;.gh altl tude satellites. The Etalon satellites offer great potential 

. 9h altitUde and vislbilil * to a larger 
stations at one time. However, the Etalon satellites have had a low 
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global ranging priority so they are rarely ranged by more than two 
stations at any one time. An intense Etalon campaign was performed 
during May and June of 1992 to determine how well these satellites can 
be used for Crustal Dynamics research and geometric data analysis. 
Unfortunately, because of the low Etalon priorities, the campaign 
yielded no simultaneous ranging data sets from 4 stations and only 1 
data set where 3 station ranged simultaneously. LAGEOS 2, scheduled for 
launch in October 1992, offers the next best opportunity to obtain a 
large simultaneous ranging data set. In addition, at an inclination of 
about 65 degrees, there will be an opportunity to obtain simultaneous 
tranB-atlantic data sets on a high priority SLR satellite. 

In the future, we expect that simultaneous geometric analysis 
of data from LAGEOS, LAGEOS 2 , and Etalonl and Etalon 2, spread over a 
period of a few days, will accurately determine baselines and 
velocities. Future modifications to the geometric analysis software are 
expected to include the additional determination of range and time 
biases from individual data sets . 
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JUNE 1988 LflGEOS POINTS 
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148°5 W 133 5 'W 118 5 'H 103 5 'M 
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